Compact Waves in Microscopic Nonlinear Diffusion 
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We analyze the spread of a localized peak of energy into vacuum for nonlinear diffusive pro- 
cesses. In contrast with standard diffusion, the nonlinearity results in a compact wave with a sharp 
front separating the perturbed region from vacuum. In d spatial dimensions, the front advances as 
ji/{2+tia) according to hydrodynamics, with a the nonlinearity exponent. We show that fluctuations 
in the front position grow as ~ f^rj, where fj, < 2+da ^ '^^^ exponent that we measure and is a 
random variable whose distribution we characterize. Fluctuating corrections to hydrodynamic pro- 
files give rise to an excess penetration into vacuum, revealing scaling behaviors and robust features. 
We also examine the discharge of a nonlinear rarefaction wave into vacuum. Our results suggest the 
existence of universal scaling behaviors at the fluctuating level in nonlinear diffusion. 
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Random walks, Brownian motions, and their contin- 
uum descriptions in terms of diffusion equations underlie 
uncountable natural phenomena [Ij [2] . Our experience 
with diffusion strongly influences our intuition. For in- 
stance, a peculiar property of the diffusion equation, viz. 
the instantaneous propagation of perturbations, appears 
puzzling when we first learn it, but eventually we regard 
it as a general property of parabolic partial differential 
equations (PDEs) which distinguishes them from hyper- 
bolic PDEs. It then comes as a surprise that the instan- 
taneous propagation of perturbations is lost for certain 
parabolic PDEs, so that phenomenologically they start to 
resemble hyperbolic PDEs. A simple class of parabolic 
PDEs which exhibits such a property is 

dtp^v ■{p'-vp) (1) 

For such equations, the spread of a initially localized pro- 
file into vacuum is not instantaneous whenever a > 0. 
More precisely, in d dimensions the front advances as 

~ ^i/(2-i-!ia). ^ ^ j,^^ |.]-^g medium is still in the vac- 
uum state [3H5] . There is no contradiction with standard 
lore as Eq. Q is non-linear. Further, this spectacular ef- 
fect occurs only when we consider the spread into vacuum 
(where diffusion coefficient = vanishes). 

Nonlinear parabolic PDEs form a fertile research area. 
We emphasize that nonlinear parabolic PDEs similar to 
Eq. ([I]) and their microscopic brethren are by no means 
pathological. For instance, Navier-Stokes equations are 
nonlinear parabolic PDEs, and e.g. in heat conduction 
the coefficient of thermal conductivity always depends 
on temperature [SI |4j, for instance it scales as ^/T for 
the hard-sphere gas indicating that the mathematical de- 
scription is similar to Eq. ([I]) with a = ^. 

Here we investigate microscopic stochastic processes 
which admit a macroscopic description in terms of 
Eq. ([T]). Amending microscopic rules, it is usually possi- 
ble to vary the exponent a. For concreteness, we focus on 
a specific exponent a — 1. More importantly, we consider 
only simple microscopic models as our major goal is to 
go beyond hydrodynamics and to analyze fluctuations. 



Two microscopic processes that meet above require- 
ments are the Averaging Process (AP) and the Ran- 
dom Exchange Process (REP). These are continuous- 
time Markov models that can be defined on an arbitrary 
lattice. At a coarse-grained level these two models mimic 
the physics of systems which are characterized by a single 
locally-conserved field. By definition, each site i carries a 
scalar variable pi > which we call energy. For concrete- 
ness, we start with the one-dimensional (Id) lattice. For 
the AP, the evolution proceeds by choosing a bond at ran- 
dom, say (i, i + and averaging the energies of the sites 
connected by the bond, {pi,pi+i) — >■ (^Ej, ^E^), with 
T,i = Pi + pi+i- For the REP, the pair energies are ran- 
domly redistributed, (pi,pi+i) — >■ (pE^, (1 — p)Ei), where 
p G [Oil] is a random number chosen uniformly on the 
interval [0, 1] in each "collision" event. If collisions occur 
independently on the energies, the above stochastic pro- 
cesses are one of the simplest microscopic implementa- 
tions of the diffusion process |7] , while energy-dependent 
collision rates yield nonlinear diffusion. 

Consider for concreteness the AP. If the collision rate 
is proportional to the total pair energy E^, the energy 
Pi {t + dt) at site i after an infinitesimal time increment 
dt will become ^Ej_i with probability Ei_idt, ^Ej with 
probability Y^idt, or will remain as Pi{t) with probability 
1 — (Ei_i -I- Y,i)dt. Averaging and taking the limit dt ^ Q 
we obtain 

'^-^^\{{pU)-'2{pI) + {pU))^ (2) 

where we write in short {pi) = (Pj(^)>, (pf) = 
etc. The governing equation ^ for average energies 
involves second-order averages, whose time evolution is 
controlled in turn by higher-order averages, giving rise 
to an infinite hierarchy which is not tractable in gen- 
eral. This hierarchy can be closed however by assuming 
that higher-order averages factorize [5], i.e. {pj) = (pi)^- 
This conjecture, a central tenet of the hydrodynamic ap- 
proach, is expected to become asymptotically valid in the 
long time limit where evolution is dominated by the 
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slowly-varying, locally-conserved hydrodynamic field. In 
this situation the spatial variance is very small, so that we 
can replace the discrete (second-order) derivative on the 
right-hand side of (2]) by a continuous derivative. Chang- 
ing now the spatial variable i — >■ a; to emphasize that we 
are employing the continuum approximation, and defin- 
ing p{x,t) = {px{t)), we arrive at 



dtp = dx{pd.j,p) . 



(3) 



i.e. Eq. ([T]) with exponent a — 1. [For the REP we 
similarly obtain dtp = 2dx{pdxp).] Generally the AP 
and REP processes with collision rates proportional to a 
power of the pair energy, T,f, are the examples [H |9] of 
microscopic stochastic processes whose hydrodynamics is 
governed by the generalized nonlinear diffusion equation 
([T]). Needless to say, fluctuating behaviors are outside 
the realm of the hydrodynamic approach. 

First, we consider the evolution of a localized peak of 
energy: p{x,t = 0) = poSix), where po = dx p{x,t) 
is the total energy which remains constant. The lack of 
characteristic spatial and temporal scales in our initial- 
value problem suggests to apply the transformation x — > 
ax, t — >■ /?<, p — >■ 7/9. Equation ([s]) is invariant under this 
transformation when — fS^y, and the conservation of 
energy leads to another constraint aj — 1 . Therefore a — 



T 



(3^^^, and our problem remains invariant under the 



one-parameter transformation group x — /3 ' x, t — 
Pt, p — >■ (3~^^^p, which suggests to seek a solution in a 
self-similar form 



p{x,t) 



Po 



(4) 



Plugging ansatz Q into the governing equation and in- 
tegrating once we flnd that the scaled energy proflle <& 
satisfies <I>(<I>' + i^) = 0, where prime denotes the dif- 
ferentiation with respect to ^. This equation reduces to 
= — ^/3 in a bump region |^| < and to $ = 
outside the bump: 




1^1 > eo 



(5) 



with A = 1/6. The conservation of energy / d^^iO — 1 
fixes the boundaries of the bump: — (9/2)^/^. (For 
the REP, due to the extra factor 2 in its hydrodynamic 
equation, A — 1/12 and = 9^^^.) The most strik- 
ing feature of the (parabolic) energy profile ([5| is that 
it is compact — the energy density vanishes at a finite 
distance from the initial peak. This contrasts with stan- 
dard expectations based on the behavior of solutions to 
the linear diffusion equation where the influence prop- 
agates instantaneously — when diffusion is non-linear, 
the 'knowledge' about the peak at the origin is accessible 

only on distances |x| < {pot)^^^ ■ 

The self-similar solution Q-IS]) was derived in [TOllllj : 
similar solutions were found earlier in studies of the inter- 
mediate stage of nuclear explosions. These solutions have 
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FIG. 1: (Color online) Right front position (0)i its stan- 
dard deviation (□) and correlation between opposite fronts 
(0) measured for the AP (blue, open) and the REP (red, 
filled) models. Inset; Distribution of front fluctuations for 
the AP and the REP. 



been investigated in the realm of hydrodynamics, which 
says nothing about fluctuations. In what follows we ana- 
lyze fluctuating corrections to the average hydrodynamic 
behavior using both the AP and REP stochastic models. 
In particular, we run extensive simulations of both mod- 
els, averaging over many different realizations with fixed 
initial condition piit = 0) = poSifl. 

We first consider the right front position, represented 
at the microscopic level by the position (t) of the right- 
most site where the energy at time t is positive. The front 
position £-f-{t) is a random quantity. Conjecturally, 



i+it)^CoiPot)'^'' + {potrv 



(6) 



where < | is a new exponent and 77 is a random vari- 
able. Figure [l] shows the measured front position and 
its standard deviations a-+{t). This confirms that the 
front propagates as t^/^ for both the AP and the REP. 
Interestingly, front fluctuations are characterized by an 
exponent p = ^ for the REP model, while p = for 
the AP, i.e. much smaller fluctuations in the latter than 
in the former for long times. This can be attributed 
to the additional layer of stochasticity (the random en- 
ergy exchange) which characterizes the REP. We also 
measured the correlation between left and right fronts, 
K{t) = [{i+i-) - (f+)(£_)](t). Opposite fronts turn out 
to be anti-correlated, K(t) < 0, due to energy conserva- 
tion, see Fig.hjwhere we plot c{t) = yJ—K{t). Moreover, 
we observe c(I) ^ with the same exponents p as we 
measured above for the AP and the REP. 

The probability distribution '^{vi) which characterizes 
front fluctuations is unknown, though its asymptotic be- 
haviors can be understood heuristically. To estimate 
the T] — > —00 tail let us flrst determine the probabil- 
ity that ^+(t) = 0. Since £4. can only increase, the event 
t+{t) = implies that during the time interval (0,t), 
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logarithm of the total probability of the above process 
scales as liin5=o2~-' ^ —i^ ~ —7^2/(1-^)^ Combining 
both asymptotics, we arrive at 
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FIG. 2: (Color online) Top: Scaled energy profiles for dif- 
ferent times in the AP model and the theoretical prediction 
$(^). The inset shows the microscopic penetration into vac- 
uum. Bottom: Collapse of the scaled excess penetration into 
vacuum for different times in both the AP and the REP. Inset: 
The scaled excess penetration decays as a (3/2)-stretched ex- 
ponential in both models. Data for the REP have been shifted 
vertically for clarity. 



the pair (0, 1) has never collided. The probability that 
this pair does not collide in an infinitesimal time inter- 
val [f, t' + dt') is [1 - p{x = 0, t')dt'] « exp[-p(0, t')dt% 
so lnProb(^+(i) = 0) = — /gdrp(0,T). Since we are 
not constraining the motion of the left front, the wave 
span still grows as t^f^ and hence the energy at origin 
decays as p(0,t) ^ due to energy conservation, 

which implies that lnProb(£+(<) = 0) ^ —t^/'^. From 
^ we see that — corresponds to a fluctuation 

77 ~ — i(i-3A')/3^ This in conjunction with above estimate 
yields In 9(77) ^ — Ijyp/li-SM) when 77 — >• —00. 

To estimate J'(?7) in the opposite i] ^ 00 limit, it is 
convenient to consider a discrete-time version of the AP 
model with pQ — 1, namely we pick up a pair of neighbors 
and make a collision with probability equal to the pair 
energy. We want to transfer energy to the right as far as 
possible. To this end we select the front pair l+£+) 
at each step j — £+ > for a collision, even though the 
pair collision rate is decaying as 2~^ , i.e. as the energy 
at the rightmost site decreases. The front moves ballisti- 
cally, £+{t) = t, which implies r] ~ t^~i^ via (|6]), and the 
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These estimates rely on the validity of Eq. ([6| for large 
values of the random fluctuation where ry diverges as a 
power of time. This assumption is known to be valid in 
numerous examples [5], where matching to the extreme 
typically allows to determine the tail scaling. 

We measured ^(ry) for the AP and REP models, see 
inset to Fig. [ij Data are consistent with a faster decay 
of negative front fluctuations for the REP as compared 
to the AP (In ^(77) ^ —\r]\'^ vs — I77P as 77 — ?> —00, respec- 
tively), although our extensive simulations did not allow 
us to access the rare-event tails where the asymptotic 
behavior ([7| becomes clear. The probability of positive 
fluctuations is closer for the REP and the AP, as ex- 
pected. Novel Monte Carlo methods allowing to probe 
rare events may help in measuring these tails [HI [T3] . 

Top panel in Fig. [2] shows the scaled energy profile 
measured for different times in AP model, together with 
$(0 of Eq. ^. The agreement is excellent for inter- 
mediate and long times {t > 10^). Moreover, fluctuating 
corrections to the hydrodynamic proflle can be measured, 
see top inset in Fig. [2] revealing a microscopic penetra- 
tion into vacuum of energy profiles which exhibits scaling 
behavior and robust features. We define the excess pen- 
etration into vacuum as A$ = t^/^{p{x,t)) — ^{xt~^/^). 
Data for different times nicely collapse onto a single curve 
under the scaling A$(^) = i-(i-3A')/3g:-[(^_^o)t(i-3M)/3]^ 

see bottom panel in Fig. [21 with 3^(z) a scaling function 
slightly different for the AP and REP models. Despite 
this difference, bottom inset to Fig. [2] shows that excess 
penetration into vacuum decays as a (3/2)-stretched ex- 
ponential in both models, i.e. 3^{z) ~ exp[— (z/x)'^/^] for 
z > 0, where x is a scaled penetration length. 

The generalization of the AP and the REP models to 
higher dimensions is straightforward. In this case the 
governing hydrodynamic equation is equivalent to ([T]) 
with a = 1. A localized energy pulse evolves again into a 
self-similar wave with compact support characterized by 
a parabolic scaled profile similar to (|5| , see Fig. |3j whose 
front advances as r/(t) = ^o(po^)^'^ + [poiYv- Front 
fluctuations in 2d are strongly suppressed: p, w 1/30 
for the REP, and again /i = for the AP. In 2d we 
also observe a similar excess penetration into vacuum of 
energy profiles, characterized by a scaling function with 
(3/ 2) -stretched exponential decay as in Id, and a scaled 
penetration length x which decreases with d. This obser- 
vation, together with the small value of exponent p, sug- 
gests that fluctuating corrections to hydrodynamic be- 
havior become less pronounced in high dimensions. 

We return now to Id and consider the discharge of a 
step initial profile into vacuum |14j . Such initial con- 
dition evolves into a self-similar, compact rarefaction 
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FIG. 3: (Color online) Space-time evolution of the scaled en- 
ergy for the 2d AP model. Inset: Convergence of the scaled 
radial profile for different times to the hydrodynamic solution. 
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FIG. 4: (Color online) Scaled rarefaction wave profile for the 
AP (red) and the REP (blue) models. Top inset: Asymptotic 
profile behavior for ^ — )■ — cxd in the AP. Bottom inset: Front 
position (O) and its standard deviation (□). 

wave. In the hydrodynamic approximation we need to 
solve Eq. ^ subject to p{x, i = 0) = 1 for a; < 
and p{x, t = 0) — for x > 0. The transformation 
X — ^/]3 X, t ^ (3t, p ^ p leaves the initial- value prob- 
lem invariant. The relevant self-similar solution is there- 
fore p{x,t) = ^(^) with = x/y/i, and the scaling func- 
tion obeys vl/xj/" -|- (vl/')2 _|_ i^vj/' = q. This equation does 
not admit an analytical solution, though it is still possi- 
ble to understand the qualitative behavior of The 



initial condition implies that ^ — )• 1 when ^ — — cxd and 
\1/ — >■ when ^ — >■ oo. An assumption that the solution 
spreads into the vacuum ad infinitum leads to a contra- 
diction '3]. Hence similarly to the compact solution ([5|, 
the energy profile vanishes at a finite value of the scaling 
variable £, = £,o, see Fig. |4j An asymptotic analysis in 
the proximity of the front shows that ~ ^Co('Co ~ 
as ^ — ^0- In the opposite ^ — — oo limit we employ a 
perturbative analysis to yield 1 — ~ e'^'^'^ [5]. This 
asymptotic behavior, confirmed in the top inset to Fig. [4j 
is not surprising since in the far-left region the governing 
equation is essentially a linear diffusion equation. De- 
spite the overall diffusive scaling, there is a well-defined 
front separating the perturbed region from the vacuum, 
which propagates in time as i+{t) = ^oVt + t^V with rj 
again a random variable. Bottom inset in Fig. |4] shows 
that (5 = 1/4 for the REP and 5 = for the AP. We can 
analyze in equivalent terms the tails of the distribution 
y{ri) of front fluctuations [9l[T4]. Furthermore, the scaled 
rarefaction waves in Fig. [4] also exhibit microscopic pen- 
etration into vacuum, characterized by scaling functions 
J'(z) with (3/2)-stretched exponential tail as in Fig. [2] 

Our results suggest the existence of robust scaling be- 
haviors at the fluctuating level in nonlinear diffusion 
models. This is remarkable since, despite the macro- 
scopic similarity between the AP and REP models, there 
is a crucial difference — the REP is time- reversible (i.e., 
it obeys the detailed balance condition), while the AP 
lacks microscopic reversibility as collisions cannot be re- 
versed. The observed signs of universality are thus in- 
triguing and deserve further study. The class of micro- 
scopic diffusive models introduced here represent an ideal 
benchmark to investigate strongly nonlinear phenomena 
far from equilibrium, and can be generalized to describe 
nonlinear driven dissipative media [8]. The simplicity 
of these models allows to go beyond hydrodynamics to 
study fluctuations in nonlinear, non-equilibrium settings. 
This is particularly interesting at the light of the recent 
breakthroughs in non-equilibrium physics within the Hy- 
drodynamic Fluctuation Theory formalism |131 115j . In 
fact, it would be interesting to study rare fluctuations 
and large-deviation properties for the nonlinear diffusive 
models introduced here, in particular in a dynamical set- 
ting as the compact wave propagation problem 0, thus 
complementing recent results on front and current fluc- 
tuations for shocks in simple exclusion processes [T4| . 
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